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Abstract. We consider an approach to the hyperboloidal evolution problem based on the 
Einstein equations written for a rescaled metric. It is shown that a confomial scale factor can 
be freely prescribed a priori in terms of coordinates in a well-posed hyperboloidal initial value 
problem such that the location of null infinity is independent of the time coordinate. With an 
appropriate choice of a single gauge source function each of the formally singulai' conformal 
source terms in the equations attains a regular hmit at null infinity. The suggested approach 
could be beneficial in numerical relativity for both wave extraction and outer boundary 
treatment. 



PACS numbers: 04.20.Ha, 04.25.D-,04.20.Ex, 04.30.-w 

1. Introduction 

An important problem in general relativity is the calculation of gravitational radiation emitted 
by self-gravitating astrophysical sources. In the isolated system idealization of such objects, 
one attaches to the far-field zone an asymptotic region in which the spacetime becomes 
flat in a certain sense. The resulting models are commonly referred to as asymptotically 
flat spacetimes. Because gravitational energy is not localizable and there is no generally 
satisfactory definition of quasi-local energy available fTl, the concept of gravitational radiation 
is rigorously defined in these models only at null infinity |2 3 1. As a consequence, one needs 
global access to the spacetime solution to discuss gravitational radiation in an unambiguous 
way. In numerical calculations, however, one typically truncates the solution domain by 
introducing an artificial timelike outer boundary into the spacetime. This practice introduces 
certain well-known conceptual and practical difficulties ||4l|5]|6l|7][8l|9l- 

A clean way to calculate gravitational radiation is to include null infinity in the 
computational domain. This can be achieved under certain conditions by employing the 
conformal compactification technique introduced by Penrose which allows the study of null 
infinity by local differential geometry lT0lfmfT2L This framework has been used extensively 
in the mathematical discussion of isolated systems ifTsl [T4l [Tsl [161 [TtI . In numerical 
applications, it has been successfully applied within the characteristic approach to spacetimes 
which can be foliated by null hypersurfaces in a regular way ifTSl [T9l l20l 1211 122]| . In regions 
of highly dynamical, strong gravitational fields, however, characteristic foliations are not well 
behaved due to formation of caustics in bundles of light rays generating the null hypersurfaces 
||23]| . While there are promising suggestions on how to deal with this difficulty such as the 
Cauchy-characteristic matching 1 24 25 , 26 , 27 1 , a successful implementation has not yet been 
achieved. 

A more general method to include null infinity in the solution domain has been suggested 
by Friedrich ll28l . In this approach, one solves Cauchy problems based on spacelike foHations 
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that intersect null infinity. Such foliations are called hyperboloidal as their asymptotic 
behavior is similar to that of the standard hyperboloids in Minkowski spacetime. While 
one solves a Cauchy problem in the sense of partial differential equations, the underlying 
surfaces are not Cauchy surfaces in the sense of differential geometry because their domain 
of dependence does not cover the entire spacetime. To avoid misunderstanding, the related 
problem is called the hyperboloidal problem. The hyperboloidal approach is promising as the 
underlying surfaces combine favorable properties of standard and characteristic approaches. 
They are as flexible as Cauchy -type foliations commonly used in numerical relativity and they 
approach null infinity thus enabling a clean treatment of gravitational radiation. 

In 1281 . the Einstein equations are reformulated for a conformally rescaled metric 
such that the equations are manifestly regular at null infinity. The resulting system of 
regular conformal field equations is larger than the Einstein equations and involves evolution 
equations for the conformal factor This approach has been very efficient in the effort to 
understand the global structure of spacetimes ifTTl |29l l30l [3TI . Numerical work on the 
hyperboloidal problem for this system has been performed mainly by Hiibner ll32l [33l [34l [35l 
|36J , Frauendiener 137l|38l|39l|40l|41J and Husa [42j|43J (see [441 for a review). Unfortunately, 
the regular conformal field equations could not yet be used numerically in the study of highly 
dynamical systems. 

In view of the recent success of numerical codes in solving the Einstein equations for a 
large variety of astrophysically interesting systems Il45ll451l47ll48ll49ll50ll . it seems desirable 
to implement the hyperboloidal problem directly for the conformal Einstein equations 
Il5ni52ll53ll54l . The idea is the following: Under a conformal rescaling of the metric g = il^g 
with a function O > 0, the Ricci tensor transforms as 

R^.A9] = R^^A9] - ^ (2 V^V,17 + Dng^,) + -^i^x^) V^f^g^.. (1) 

Here, is the Levi-Civita connection of the rescaled metric g and □ := g^'^V^Vi/. The 
Einstein vacuum field equations, G^^ [g] :— R^i, [g] — \gfj.i^R[g\ = 0, are thus equivalent to a 
similar system for the conformally rescaled metric 

Gf,M = T^,m (V^V.r! - U^g^,,) - ^{V g,,,. (2) 

To include null infinity in the computational domain, one allows the conformal factor Vt to 
vanish in a suitable way and solves the above system directly. It has been pointed out f55l 
that there are two major difficulties in this program, even if one assumes that the conformal 
extension is regular. First there is the question of how to fix f2. The conformal factor is related 
to the asymptotic structure of the spacetime, therefore it must be determined jointly with the 
metric. Secondly, there are terms involving divisions by powers of which are formally 
singular at {il = 0}. The question then is how to make sure that these terms attain regular 
limits at null infinity. 

Assuming a regular conformal extension, these difficulties have been solved in the 
characteristic case by Tamburino and Winicour |56|. Their solution includes a certain choice 
of coordinate and conformal gauge in which the formally singular terms in (|2]i attain regular 
limits at null infinity. Their construction provides, in a sense, the mathematical basis for 
characteristic codes that include null infinity in the computational domain. The problem in 
the hyperboloidal case has been open. 

We address the above mentioned difficulties in two propositions. The first one treats 
the case of a strictly positive conformal factor and states that it can be prescribed almost 
arbitrarily in terms of local spacetime coordinates in a well-posed Cauchy problem for the 
Einstein equations. For the second proposition we allow the conformal factor to vanish in a 
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certain way so that null infinity can be included in the computational domain. We see that a 
class of gauges can be chosen in a hyperboloidal initial value problem in such a way that the 
formally singular terms due to conformal compactification attain regular limits at null infinity 
for a certain class of asymptotically flat spacetimes. 

2. A hyperbolic reduction for a positive conformal factor 

The system (|2|i has the form of Einstein equations with source terms. In general, such 
a system must be completed by additional equations derived from the Bianchi identities, 
V^G^^ = VT^i/ = 0, implying equations of motion for the source functions. In our 
case, however, there are no additional equations required for the conformal factor VL as shown 
below. 

Proposition 1; The conformal Einstein equations (O admit a well-posed initial value 
problem for an arbitrary Q. £ C^(K'^, ]R_|_). 

First we show that the Bianchi identities are satisfied for a positive, sufficiently 
differentiable conformal factor. We calculate 

V^r^,[r!]= -j^V'^r! (2V^V,f] + m.g^,)- 

- ^ - v,no) + ^ x^)^^^- (3) 

Contracting the commutation relation V xV ^ Vpil — Vi/V a Vpil — R^^j^Vcr^ with g^P and 
exchanging derivatives we get DV^ri — ViyOr? — RJWail. Using this relation, equation Q 
with the definition of the Einstein tensor Gfj,^[g] and the conformal source tensor T^j^[il], we 
get 

V'Tp.p =-^V^n (G^,[5] - Tp.P) = 0. (4) 

We see that the Bianchi identities are satisfied by virtue of the conformal Einstein equations 
for a non-vanishing ft that is at least three times differentiable, i.e. ft E C^(M^, but is 
otherwise arbitrary. Therefore, il can be regarded as a free function. We can write down some 
suitable equation for it consistent with the above calculation or prescribe it directly in terms 
of some yet unspecified coordinate system as long as f2 7^ 0. 

The next step is to set up a well-posed initial value problem for (|2|i. We employ 
the general wave gauge reduction of the Einstein equations, also known as the generalized 
harmonic reduction for historical reasons, where the argument of well-posedness for ^ can 
be taken almost directly from ll57ll58l with only a minor modification. 

Regarding the Ricci tensor R^'^ as a differential operator acting on the metric g, we can 
write the conformal Einstein equations in a local coordinate system {a;^}p=o.i,2,3 as 

R""" [9] = Ig^'dxdpg'^" + v^^r'') - g^^g^^ r^^^r;, = t^-^ ~ ig^.r, (5) 

where we have defined the contracted Christoffel symbols := g'^^F^^ = — Dx^, set 
yppz. ^ gt'P{dpT'' + V'p^T^) and T = g^^Txp. The principal part of the operator i?^"^ is 
of no known type. It was recognized by Choquet-Bruhat [59) that one can always choose a 
wave gauge, historically referred to as harmonic gauge, at least locally, so that the contracted 
Christoffel symbols vanish, F^ — —OgX^ = 0, and the system (|5]l, reduces to a quasi-linear 
system of wave equations. This reduction technique led to the first local existence result in 
general relativity 1591 . 
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The reduction based on the wave gauge was generalized to arbitrary coordinate systems 
by Friedrich with the introduction of gauge source functions II6OII . In the general wave gauge, 
the coordinates are constructed as solutions to an initial value problem for the semi-linear 
system of wave equations DgX^ = —T^^ = —F^^ with prescribed functions F'^{x,g) that 
can depend on the coordinates and the metric. These functions act as source functions for the 
coordinate gauge, hence the name gauge source functions. Note that the general wave gauge, 
in contrast to the wave gauge described above, is not a specific choice of coordinates but a 
particular way to prescribe general coordinates in an initial value problem. 

The reduced system for (|5]l is then obtained by replacing the contracted Christoffel 
symbols with the gauge source functions F^. The result is a quasi-linear system of wave 
equations for the metric components which can be written as 

G'^'ig] = T^'-'in] + V^^'C"^) - ^{VxC^)g'"', (6) 

where — — are called the constraint fields. We want to study the Cauchy problem 
for this system. We will only point out certain aspects that play a role in later considerations 
or that are different from the detailed discussion in |58 1. 

The Cauchy data on an initial hypersurface 5 = {a;'' = 0} consist of g'^'^ls and dog^'^ \s . 
Assume we are given on 5 a Riemannian metric hab and a symmetric tensor field Kab as a 
solution to the Einstein constraint equations where a,b = 1, 2, 3. We choose gauge source 
functions F^ (x'*' ) and four functions on S that correspond to initial data for the lapse function 
a > and the three components of the shift vector /?" . In the interior, these functions should 
be chosen such that do is timelike which implies — hab(i°'ff' > 0. We will later allow do 
to become null at the outer boundary (see the discussion leading to (fT2l i). We obtain the data 
g'^'^ls via the decomposition 

The data dog'^'^\s is determined such that C^l^ = and Kap is the second fundamental form 
on S. Standard theorems guarantee that we can find a unique solution to the Cauchy problem 
for the reduced equations (|6]l that depends continuously on the initial data. The solution spaces 
of ^ and (|2|l are equivalent if the constraint fields vanish. The Bianchi identity, V^G^'^ = 0, 
and dUi together with (|6]l imply the following subsidiary system for the constraint fields 

aC + R'iC" - ^V^f7 J^V^^C") - \{^\C^)9^''^ = 0. (8) 

Initial data for the evolution equations has been constructed such that C^j^ = 0. From the 
evolution equations evaluated on S it follows that SqC^Is — 0. The uniqueness of solutions 
to the Cauchy problem for the semi-linear, homogeneous system of wave equations for 
given in ^ then implies that the solution to the reduced system ^ satisfies — away 
from the initial surface S. Thus we have shown that there is a well-posed hyperbolic reduction 
for the conformal Einstein equations 

We shall briefly elaborate on how the free prescription of the conformal factor is to be 
understood. We can not prescribe the conformal factor as a function on the manifold because, 
in an initial value problem, we do not know the manifold. The prescription of a function for 
the conformal factor determines only its representation in terms of coordinates which are yet 
to be constructed during the solution process. Invariant properties of the resulting conformal 
factor will depend on initial data and the choice of gauge source functions. The essential 
property of (|2|i that is responsible for this feature is its conformal invariance, in the sense that 
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if (A^,g,i7) is a solution to £'^,^[5] = rp,y[ri], then {M. , uP' g ^ ujVI) with a positive function cj 
is a solution to G^iyluj^g] = T^i,[loVI]. The system (|2]i determines the conformal class of the 
metric g in contrast to the Einstein equations which determine the isometry class of g. This 
allows us to prescribe an arbitrary coordinate representation for the conformal factor as long 
as there are no geometric requirements to be satisfied. 



3. Choice of gauge at null infinity 

In this section, we deal with the formally singular terms in the conformal source tensor 
T^^[ri]. We first present a preferred conformal gauge at null infinity in which each of the 
conformal source terms attains a regular limit at null infinity in a given conformal extension. 
Then we present how this choice of gauge can be achieved by a suitable choice of gauge 
source functions in a hyperboloidal initial value problem where the conformal factor has been 
prescribed explicitly in terms of coordinates. We restrict our attention to future null infinity, 
denoted by The treatment of past null infinity follows by time reversal. 



3.1. The preferred conformal gauge at J'^ 

Assume that a solution {^M. , g) to the Einstein vacuum field equations has been given which 
admits a conformal extension 17, fi) including a smooth piece of The existence of 
a broad class of such solutions is due to ll28ll6n . It has also been shown that solutions exist 
which not only admit a smooth piece of J^"*" but a complete [62] (see also I29l l63ll64l ). 
We show that ^+ C {51 = 0, rffi ^ 0} is a shear-free null surface independently of the 
conformal gauge IITSI ISSI. Multiplying Q with ^ and evaluating it along ^+ we see that 
g^f'V p^\ j: Jr = 0. This shows together with dn\j^+ ^ that ^+ is a null surface. Now 
multiply (O with fl and take its trace-free part along J^'^ to get 



The relation above is independent of the conformal gauge because we derived it from the 
conformal transformation behavior of the Einstein tensor (|2]i. Another way to see the 
conformal invariance of (|9]l is to consider the transformation behavior of (|9]l under a further 
rescaling of the conformal metric given by 

g'=Lj'^g, n'^ujn, uj>OonM. (10) 

We have — LuWf^W^il + g^^, V'^JIVac^. The trace of this relation reads 

= ^ (cjnsi + 4v-^f7VAw) , (11) 



whence we get 



= 0. 



To see that (|9| implies shear-freeness of J^'^, we introduce in a neighborhood of 
a null vector field that satisfies — V^il. We complete to a Newman-Penrose 

complex null tetrad {I, k,m,fh) satisfying the usual relations f66|. Newman and Penrose 
introduce twelve complex functions called spin coefficients. We are interested in two of them, 
namely a :— m^^m'^V ^l^ and p :— m^^fri'^V ^l,^. As discussed in ll66l . when Z'^ is tangent to 
an affinely parametrized null geodesic, a can be interpreted as the complex shear of the null 
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geodesic congruence given by and the expansion of the congruence is characterized by p. 
We see with ^ 



In our case, the null generators of are not necessarily geodesic, that is, in general they 
do not satisfy l^V = on However, under a rescaling of given by {l')^^ — 91^ 
with a positive function 9, the spin coefficient a transforms as a' — 9a, so the vanishing of a 
is invariant under such a rescaling which we can use to make l^^ geodesic. We conclude that 
J^'^ is a shear-free surface in any conformal gauge. 

While the vanishing of a and thus the shear-freeness of ^+ is conformally invariant, 
the vanishing of the expansion of J^+, characterized by □r2| jir+, depends on the conformal 
gauge. As can be seen from (fTTl i. given a conformal extension we can always find 
a rescaling ( fTol l such that □'r2'|jif+ = by solving the ordinary differential equation 

We call the conformal gauge in which the expansion of ^+ vanishes a preferred 
conformal gauge. This gauge has been useful in mathematical studies because of its special 
properties 16511671 . It is also the gauge choice of Tamburino and Winicour [56|. By (|9]l and 
(O we see that in a preferred conformal gauge 



which implies that each conformal source term in (|2]i attains a regular limit at J^"*". 
3.2. The choice of gauge source functions at J''^ 

The above construction of a preferred conformal gauge assumes that a conformal extension 
has been given. We are interested, however, in the case where only initial data has been given 
and we would like to know how to choose the gauge source functions suitably so that the 
hyperboloidal evolution is performed in a preferred conformal gauge. 

Proposition 2; Assume hyperboloidal initial data has been given whose evolution admits 
a smooth conformal compactification at J^+. A preferred conformal gauge can be achieved in 
a general wave gauge reduction o/(|2| by using as a coordinate near {fl — 0} and choosing 
the related gauge source function such that F^\_^^2=q} — 0- 

We can use the conformal factor as a coordinate near {fi — 0} because df^|{si=o} 0- 
In the general wave gauge, there is a gauge source function related to each coordinate via 
DgX^ = —F^. Using the conformal factor as a coordinate, we see that the value of the gauge 
source function F^^ — — Dil at is a direct measure of the expansion of J^+. Setting it to 
zero makes the expansion vanish and hence gives a preferred conformal gauge. 

While its proof is very simple, proposition 2 is a remarkable property of the conformal 
Einstein equations. In general, it is not known how a desirable gauge can be achieved by a 
suitable choice of gauge source functions |68| because geometric properties of coordinates 
depend not only on this choice but also on initial data in an essentially non-linear way. The 
above proposition is very special in this respect as it states that one can, by a suitable choice 
of a single gauge source function, fix an a priori known coordinate surface to be a null surface 
free of shear and expansion that corresponds to null infinity. 

In a practical numerical calculation one may proceed as follows. One calculates on a 
three dimensional surface S with a two dimensional boundary E hyperboloidal initial data 



m^m-^ 5^,nr! = 0, p\,_ 




VpV^rj|,y+ = 0, and \im - g^'^S/ x^'s7 = 0, 
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{S,hab,Kab) whose evolution admits a smooth conformal compactification at The 
existence of such data has been studied in ||6TI . In the case of isotropic extrinsic curvature, 
that is, under the assumption Kab — "j hab, the data needs to be such that the trace-free part 
of the second fundamental form induced by hab on E vanishes. This condition corresponds 
to shear- freeness of J^+. The case with more general data has been considered in 1691 iTOl . 
One can construct the data in such a way that the conformal factor has the form ft — 1 ~ r 
in a neighborhood of E where r is a coordinate whose level surfaces have spherical topology 
II52II . The form of the conformal factor in the interior can be chosen freely to satisfy practical 
needs, for example to obtain high resolution in certain domains. Initial data for lapse and shift 
should be chosen such that = — |r=i — do — Sq, which implies with ^ 

a\^ = Vh^, (3-\j, = -h''-. (12) 

The gauge source functions can be given quite freely as long as the condition F^\r=i ~ is 
satisfied. Then each of the conformal source terms in 

2 

rp /-pr rr 

will attain a regular limit at ^+ and the system (|2]i can, in principle, be solved. In practice, 
however, the accurate calculation of these limits poses a major computational difficulty that 
needs further study. 

4. Discussion 

We studied the hyperboloidal evolution problem for the conformal Einstein equations. Two 
propositions have been made that suggest an alternative method for numerical calculations of 
asymptotically simple spacetimes including null infinity, substantiating earlier expectations in 
this direction I52l |5^ l54l . 

Proposition 1 states that one can formulate a well-posed initial value problem for the 
Einstein equations in which a positive conformal factor is prescribed freely in terms of 
coordinates. In a sense, this result can be regarded as trivial because a rescaling with a positive 
factor essentially amounts to a relabeling of the metric. On the other hand, the possibility 
of prescribing the conformal factor in a well-posed initial value problem is very convenient 
because it allows us to set up a formulation of the Einstein equations for a rescaled metric 
in which the conformal factor is independent of the time coordinate in the whole solution 
domain. Depending on practical demands, one can choose the conformal factor in certain 
regions to be unity so that the standard Einstein equations are obtained, or larger than unity 
for refinement. 

Proposition 2 shows how one can construct the solution to a hyperboloidal initial value 
problem directly in a preferred conformal gauge by a suitable choice of a single gauge source 
function. In this gauge the conformal source terms attain regular limits at null infinity. In 
combination with Proposition 1, this allows us to include null infinity in the solution domain 
in such a way that the outer boundary coincides with null infinity. It is a remarkable feature 
of the conformal Einstein equations in the general wave gauge that the constraint equations, 
the evolution system, the gauge system and the subsidiary system all work in harmony with 
conformal transformation formulae allowing us to make the geometric properties of null 
infinity manifest on an a priori known coordinate surface. 

We have assumed vacuum throughout the calculations. Matter can be included by solving 
the matter equations with respect to the conformally rescaled metric if the matter fields remain 
confined to a finite region of space or fall off sufficiently fast towards infinity. 



Hyperboloidal evolution with the Einstein equations 



8 



The approach presented in this article has various advantages that are favorable for 
numerical calculations. First of all, the method is tailored for the treatment of the asymptotic 
region such that the interior scheme is not necessarily modified. In fact, hyperboloidal 
foliations can be made to coincide with Cauchy-type foliations in the interior where the 
conformal factor can be set to unity if desired 171]. In addition, the principal part of the 
equations to be solved in the exterior domain is identical with the standard Einstein equations 
which we know how to treat numerically for a large class of astrophysically interesting 
configurations. Furthermore, the location of null infinity is known a priori simplifying 
numerical outer boundary treatment and wave extraction significantly. The calculation of the 
news function is in our coordinates very simple, because the solution is obtained directly in 
Bondi coordinates at null infinity in contrast to the characteristic approach where the freedom 
in the null coordinates is fixed in the interior of the spacetime 1.721 l73l . 

A successful implementation of the hyperboloidal approach would solve two major 
problems in numerical relativity, namely the outer boundary problem and the wave 
extraction problem. One should be aware, however, that many promising ideas appeared 
during the history of numerical relativity which encountered major difficulties in practical 
implementation that remain unsolved. In our case the most delicate issue is whether the 
formally singular terms can be numerically calculated in a stable manner. This question 
is open and is left for future work. We emphasize, however, that the viability of a similar 
calculation has been demonstrated within the characteristic approach long time ago [74|. 

It should be pointed out that there is a theoretical limitation to the hyperboloidal 
approach. It does not allow us to calculate global spacetimes because spatial infinity is 
not included in the computational domain. A detailed study of the asymptotic behavior of 
gravitational fields near spatial infinity is difficult and an active area of research Il64ll75ll76ll77l 
|28J. It seems that one should first achieve a successful implementation of the hyperboloidal 
problem for astrophysically interesting configurations before attacking the more challenging 
problem of spatial infinity. 
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